Corrections to the Newtonian potential 
in the two-brane Randall-Sundrum model 
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We calculate the Newtonian potential in the two-brane Randall-Sundrum model, emphasizing the 
effect of the finite distance between the two branes. The result obtained is quite natural: When 
the distance in the potential is small compared to the brane separation the two-brane model is 
indistinguishable from the one-brane model, whereas when the distance is large the bulk dimension 
behaves like an ordinary compact dimension, with an exponentially decreasing correction to the 
four-dimensional potential. The contribution from the radion is also included, and is found to give 
only a multiplicative factor of order 1 in the correction. 

PACS numbers: 04.50.+h, 04.62.-fv, ll.10.Kk, 98.80. Jk 



I. INTRODUCTION 

In this paper we study the Randall-Sundrum I model 
0, which consists of two 3- branes embedded in five- 
dimensional anti de Sitter space. Our main focus will 
be the static gravitational potential between two point 
particles on the visible brane, and in particular the effect 
the finite gap between the two branes has on the poten- 
tial compared to the Randall-Sundrum II model |2| with 
only a single brane. 

The two-brane model is mostly used in connection with 
particle physics, where it gives a natural explanation 
of the hierarchy between the electroweak scale and the 
Planck scale, without the need to introduce very large 
(or very small) dimensionless parameters in the theory 
1] . For a survey of the phenomenology of theories with 
extra dimensions, see e.g. Q and references therein. As 
for the problem of graviton localization and the deriva- 
tion of the gravitational potential, most of the recent 
activity has focused on the one-brane model 0, Q . How- 
ever, some work has also been done with two branes in 
this context, see e.g. |EI3' an d 111 where induced curva- 
ture terms in the actions for the two branes are consid- 
ered. The cosmology of the two-brane model is studied 
in e.g. 0, and a mechanism for stabilizing the brane sep- 
aration is discussed in [T^j. In this paper we will simply 
assume that the branes are static and fixed at a constant 
distance. 

Before doing any calculations, we should be able to 
give a qualitative description of what the potential should 
look like. The deciding factor must be the ratio between 
the distance r in the potential and the separation between 
the two branes. For short distances the second brane is 
far away, relatively speaking, and should therefore be in- 
visible to the potential. The second brane could just as 
well be sent away to infinity, meaning that the potential 
must be the same as in the one-brane model for short dis- 



tances. Of course, if the distance r is very short we can 
also ignore the entire five-dimensional curvature, mean- 
ing that V(r) ~ 1/r 2 when r — > 0. On the other hand, for 
large distances the more detailed structure of the fifth di- 
mension must be invisible to the potential, and we should 
get the same result as with an ordinary, compact extra 
dimension. That is, the potential should essentially be 
four-dimensional, V(r) ~ 1/r, with an exponentially de- 
creasing correction from the fifth dimension. 

However, it is not the physical distance y r to the second 
brane (as measured with standard rods) that sets this 
scale, but rather the conformal distance z r ~ e Ky, \ This 
can be understood from the fact that all particle masses 
as measured on the visible brane are reduced by a factor 
g-M2/r as compared to the masses on the hidden brane 
The exponential factor originates from the warped 
geometry of the fifth dimension. Since the gravitational 
force is mediated by massive gravitons when observed 
in four dimensions, and massive propagators result in 
a factor e~ mr , the distance scale is determined by r ~ 
mT 1 ~ e+Mjfr ^ Zr ^ j e the conformal rather than the 
physical brane separation. 

This paper is organized as follows: First, the model 
is presented in section [HJ In section IIIII we summarize 
the governing equations for graviton propagation in this 
particular background, using the results of @|- We then 
proceed to find the spectrum of graviton masses in section 
IIVI allowing us to finally derive approximate expressions 
for the gravitational potential in the limit of short and 
large distances, respectively. 



II. THE MODEL 

We assume that the five-dimensional spacetime can be 
parametrized by the metric 



A 2 (y)r) f , u dx t "dx u - dy 2 , 



(1) 
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where A(y) is some function of the coordinate y of the 
fifth dimension, called the warp factor. The signature 



2 



we use is rj^v = (+1,-1,— 1,-1). The two branes are 
located at y — and y = y r with an orbifold symmetry 
S 1 jZi in the y-direction, and we assume that the only 
source of gravity is a cosmological constant A^ in the 
bulk, and a tension A and X r on the two branes. The warp 
factor A(y) is determined from the Einstein equation just 
as in the one-brane case. With the fine-tuning Kb = 
— M _6 A 2 /6, M being the five-dimensional Planck mass, 
the solution is 



A(y) = 



(2) 



where fi = y/\As\/6 = M _3 A/6. The bulk space is anti 
de Sitter since Ab < 0, and the fine-tuning means that 
the effective four-dimensional cosmological constant van- 
ishes. This configuration is often called a critical brane. 
The two branes impose the boundary conditions 



A 



■\M- 



'■x. 



A 



T M~ 3 X 



(3) 



V=Vr 



on the warp factor, where [A'] y=0 = A'(0 + ) - A'(0~) = 
2A'(0+) and [A%= Vr = A' (y+) - A' \y~) = -2A'(y-) are 
the jump discontinutics of A'{y) across the two branes. 
The boundary condition at y = simply fixes the so- 
lution whereas the condition at y = y r imposes a 
constraint on the tension A r of the second brane. With 
A(y) = for \y\ < y r we get A r = -6/xM 3 = -A, 

which is the well-known result that the two branes have 
tensions with the same magnitude but opposite sign. The 
physical brane where we are supposed to be living, is usu- 
ally taken to be the one with a positive tension, i.e. A > 
and A r < 0. 

A preliminary expression for the four-dimensional 
Planck mass, ignoring the contribution from the ra- 
dion, can be found by integrating over the y-dimension 
in the five-dimensional gravitational action S = 
— 7}M 3 f d A x dyy/gR. All five-dimensional quantities are 
here denoted with a hat. The five-dimensional curvature 
scalar R is related to the four-dimensional one through 
R = (l/A 2 )R+ . . . (replacing for a moment the flat four- 
dimensional metric r]^ v by g^v)- Using y/g = A\fg and 
demanding that S = — |Mpj J d 4 Xy/gR+ . . ., we then get 



M|i = M 3 



A 2 (y)dy 



M 3 
A 4 



f - e 



-2» 



(4) 



The four-dimensional Planck mass thus depends only 
weakly on y r as discussed in and the parameters M 
and fi can be chosen to be of the same order of magni- 
tude as Mpi. The presence of the radion changes this 
expression slightly, see eq. PU)|. 



III. GRAVITON PROPAGATION 

We study the propagation of gravitons in the model by 
looking at a perturbation 



ds 2 



A 2 (y) (ritu, + V) dx^dx" - dy 2 



(5) 



to the metric [l2|- This analysis is exactly the same for 
two branes as for a single brane, only that the regulator 
brane in p| now remains at a finite distance. We will 
therefore only summarize the main results here, referring 
to |5( for the more detailed calculations. 

Imposing the usual gauge conditions d^h^ = and 
Tl^h^v — 0, the wave equation for reduces to 

1 A A' \ 

? V 2 a -^- T 9JV = 0, (6) 

where V„ is the four-dimensional Laplacian, V„ = d a d a 
in flat space. The separation of variables h fIlV {x,y) — 
Gp V (x)${y) yields 



(Vl + m 2 ) G^{x) =0, 



<£"(y) 



AA' 
~A 



*'(!/) 



A 2 



$(y) = 0. 



(7) 
(8) 



The tensor G liV {x) describes massive spin-2 particles in 
four dimensions, and the allowed values of the mass 
m are determined from JHJ). Changing to the confor- 
mal coordinate z where dy/dz — A(y), meaning that 
I + (i\z\ — e 
reduced to 



and writing $(y) = A 3 / 2 u(z), |JSJ is 



[-d 2 z + V{z)] u(z) = m 2 u{z) 



(9) 



where 

V{z) = \{A') 2 + \AA" 



A(l + n\ Z \)> 



— 3/j,S(z) 



3//. 



f 



[iz r 



-S(z-Z r ). (10) 



(A' still means the derivative of A with respect to y.) The 
two delta functions in the Schrodinger potential result in 
the boundary conditions 



2u'(z r ) + 



2w'(0) + 3/iu(0) = 0, 
3// 



1 



flZ r 



-u{z r ) = . 



(11) 

(12) 



The solutions to 10 with the first boundary condition 
(|ll|l imposed are 

u (z) = iv (i+ M Mr 3/2 , 



,(z) = N m y/l + n\z\ 



x r 2 [^(i + M |z| 



Fi(-) 

- 7 -^j 2 [^(i + M |-l 



(13) 

where No and N m are normalization constants, and J n (x) 
and Y n (x) are Bessel functions of order n. The second 
boundary condition H12[l determines the allowed mass 
eigenvalues, and can be simplified to 



Jl(j) " Jl[^(l+^)] 



(14) 



for the massive modes. The massless mode will always 
be present. 
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IV. CORRECTIONS TO NEWTON'S LAW 



A. Finding the mass eigenvalues 



As shown in |5J, the gravitational potential between 
two point particles on the physical brane at y = 0, due 
to the exchange of virtual gravitons, can be written 



1 



rjif-Lts TjV" L J rp 

V^i r m,2 J 1 r """ si ' 



nk) =II 3^Km,0)| 2 



fJ,V p("l) rpO.fi 



k 2 



(15) 



The energy momentum tensors are Ty V = miu^u 1 ' = 

mi^g^Q and — m 2 <5Q<5g , and Pj^g is the polariza- 
tion tensor of a spin-2 particle in four dimensions with 
mass m and momentum k. The normalization constants 
in (|13[1 are found by requiring J_ r Ui(z)uj(z)dz — <5y, 



which means that 



dz 



(i + mM) 3 



1 



,—2ny r 



(16) 



Using P^ Q0 



— h and Pi™n°' = §, together with M P1 2 = 



2 

2 - 1 0000 — 3 
8irG where G is the Newtonian gravitational constant 

we get after Fourier inverting 1)15(1 



V(r) 



Gm\m2 



4(1 



■if I 



m>0 



Using the identity J n (x)Y n+ i (x) - Y n (x) J n+ i (x) = 
the expression for |u(m, 0)| 2 can be simplified to 



(17) 



2 

7TX ' 



|u(m,0)| 5 



2/i 



3 r 



m(l+:r) 



(i 



]-n(^)J 2 [ 



1 2 

|j dx 



(18) 



In order to proceed further we must find the allowed mass 
eigenvalues by solving (|T4Tl . at least approximately. 

If we also include the radion, which has been ignored 
until now, the potential changes to 



V(r) 



Gm\m,2 
r 

4 1 



-2ny r 



1 



3/i 1 



l e -2uy r 
3 m>0 



2 — mr 



So far, the conformal distance z r to the second brane 
has been considered a free parameter of arbitrary magni- 
tude. However, the original idea of Randall and Sundrum 
PJ was that the two-brane model should explain the hier- 
archy between the electroweak scale Mew ~ 1 TeV and 
the Planck scale Mpi ~ 10 18 GeV, which is obtained by 
setting 



= 1 + jlZ r 



10 



15 



(21) 



In the following, we will therefore use as a simplifying 
assumption that jiz r is a very large number. The general 
features of the results we will obtain in the end, however, 
are independent of this assumption. In appendix lAl the 
other extreme case fiz r <C 1 is studied, showing that the 
only thing that affects whether the two-brane model can 
be distinguished from the one-brane model is the ratio 
r/z r . 

The function f(x) = Y\(x) / J\{x) is almost periodic, 
and approaches tan(x — 3 j : ) fast for large x. From 114|) 
we therefore get m n z r ss nix as a first approximation. 
However, we can find a much more accurate expression 
than this for the lowest eigenvalues in the limit [iz r ^> 1. 



Since ■ 
ments: 



fJ,Z r 



« 1 we expand f(x) for small argu- 



f(f) 



O(lnS). 



(22) 



The pole at x = must coincide with a pole at x = c„, 
that is, the n-th zero of the Bessel function, Ji(c„) = 0. 
With the expansion f{c n + e n ) = Y 1 (c n )/e n J[(c n ) + 0(l) 
and demanding that /(^) = /[^(l + fj,z r )] = f(c n + e„) 
we then get 



7rm 2 Yi(c„) 



+ 0(=fln^), (23) 



4/i 2 J{(c„) " ' "' " 



and hence 

YR"n C n -T 6 n 

ft 1 + flZ r 

c n TrYi(c n )cl 



l + HZr 4J{(c„)(l + HZ r ) 3 



o 



111- 



(l+/iZ r ) 5 1+nz. 



(24) 



For large values of n we use Ji(c„) ~ cos(c„ — ^) = 0, 



(19) which means that c„ ~ (n + j)tt. 



and the four-dimensional Planck mass is then 
M 3 l-e- 2 » 1 



pi 



/i 1 + \e- 2 Mr 



8ttG 



(20) 



Comparing eq. {TH|) and with {Hi and |@}, respec- 
tively, we see that the whole effect of the radion is the 



extra factor (1 



'2^y, 



)• 



B. Short distances, fir <C 1 

For very short distances we should expect to get the 
same result as for flat, five-dimensional Minkowski space, 
i.e. we should expect the relative correction A(r) to the 
Newtonian potential to be proportional to 1/r, where 
V(r) = V (r)(l + A), V (r) = -Gmym^jr. The full 



4 



potential will then behave like V(r) ~ 1/r 2 when r — > 0, 
since A ^> 1 in this limit. 

From (|19[) we see that r — > corresponds to letting 
?ti — > oo. More precisely, from the assumption /ir -C 1 
we get m ~ r _1 3> /i, allowing us to use the asymptotic 
expressions for the Bessel functions in Ijl8(l . The integral 
thus reduces to 



ir^m Jo F 7r z TO z 



(25) 



which means that |u(m, 0)| 2 ~ l/z r . By also using m n 
rm I z r we then get 



-1jj,y r 



oo 

\ * g—nnr/zr 



3/iz r 1 + i e - 2 » ^ 

4 1 - e" 2 ^ 1 
3jii2 r 1 + \e~ 2 ^ e^/^ 



1 



(26) 



Since fir <C 1 implies r <^ z r we can expand the expo- 
nential in powers of obtaining the final result 



A 



1 



,-2^y r 



37T/Ltr 1 



-2fiy r 



fir < 1 



(27) 



The only dependence on the brane separation y r lies in 



the factor (1 — e 



-2/J.Vr 



)/(l 



), which is exactly 



the factor that enters into the five-dimensional Planck 
mass M from l|20|) . To be more spesific, since A ^> 1 the 
complete potential is 



V(r) ~ - 



4G 1 - e" 2 ^ 



1 



mim2 

67r 2 M 3 r 2' 



(28) 



by using **""")) and G = (87rA/|j) _1 . The natural way of 
defining the gravitational constant Gd in D spacetime 
dimensions is by requiring the gravitational force to be 
F = G 017111^2 / r D ~ 2 . Since the force is given by the 
derivative of the potential we thus get G5 = (37r 2 M 3 ) _1 , 
and hence V{r) ~ —Gsmimi/Zr*. 

From l (2*7)l we also notice, of course, that we reproduce 
the one-brane result by taking the limit y r —> 00. 



C. Large distances, fir ^> 1 

The limit of very large distances, r — > 00, corresponds 
to m — * (or more precisely m <SC /z). Things will there- 
fore be a little more complicated, since we can't use the 
asymptotic expressions for the Bessel functions in i|18|) 
anymore. Instead, the integral must be evaluated more 
or less analytically as it stands. Changing to the variable 



L (l + x) we get 



J^Y^-Y^Uu) 



du 



Jl{f)h - 2J 1 (*2)y 1 (s)/ 2 + Y 2 (f)I 3 . (29) 



The two integrals 1% = J uY 2 2 (u)du and h = J wJf(u)du 
are easily expressed in terms of other Bessel functions, 
whereas the integral I2 — J uY2(u)J2(u)du is a little 
more complicated and involves the Meijer G-function 
fl3| . However, all the different terms can be expanded 
in powers of mj fi, and they turn out to be of completely 
different orders of magnitude. The dominant term can 
be shown to be the upper limit u = c n in ijj, where 



2 and thus 



4/^ 2A 2 J 2 (c») m) 

7 ™ - ~2 J 1 3\ C n) - h—A ) (3U) 

7T 7TI 7T 771 

where we have used h(u) — \u 2 [Jf (u) — Jx(u)Js(u)] 
together with Ji(c„) = 0. Inserting this into l|18(l gives 



\u(m,0)f 



/JC 2 jf(c„) (1 + [iZ r ) 2 Jl(c n ) 



(31) 



For very large distances, r > z r , we get the dominant 
contribution from the first eigenvalue mi, 



A 



1 - e - 2m - 



3(1 + /iz r ) 2 J 2 (ci) l + Ie- 2 ^ 
8.21953739 



(l + /JZr) 2 



-3.83170597 r/z r 



r » z r , 



(32) 



where we have ignored the terms e~ 2 ^ Vr (since this repre- 
sents a higher order correction) and inserted the value of 
ci in the last line. The correction thus decreases exponen- 
tially for large distances, which is the familiar result for 
a compact extra dimension. The correction is, however, 
strongly suppressed compared to an ordinary compact 
dimension, because of the factor (1 + fiz r )~ 2 ~ 10 -30 . 

For intermediate distances we must sum over all the 
allowed eigenvalues. Using the approximations J|(c„) ~ 
-r 2 — and c„ ~ (n + 4)tt we then get 



2rr 2 



OO 



3(l + MZr) 2^ 

7r 2 e - 7 rr/4^ 5 _ e -nr/z 



24(l + /jz r ) 2 sinh 2 ( 



(fl+ l) e -(»+l/4W^ 

fir > 1 . (33) 



2z, 



Since we have assumed fiz r 3> 1 this expression is also 
valid for r <C z r , as long as we still have \ir 3> 1. In this 
limit (|33|l reduces to 



A 



3/i 2 r 2 ' 



1 <C [ir <C /iz r • 



(34) 



which is the same result as for a single brane and in 
agreement with |8j. 

In figure []] the two expressions *""7)l and "*"*"|) for short 
and large distances are combined into a single plot (us- 
ing a simple interpolation for the transition at /zr ~ 1). 
As is evident from the figure, there are three different 
regions, defined by the scale of the two parameters \i and 
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z r in the theory: For short distances (fir < 1) the poten- 
tial is essentially five-dimensional, just as expected, since 
spacetime can always be considered locally flat. For inter- 
mediate distances (p^ 1 < r < z r ) we have the dominant 
term A ~ 2/3/i 2 r 2 just as in the one-brane model. Here 
the distance is large enough to be affected by the warped 
geometry of the fifth dimension, but still too short to no- 
tice the presence of the second brane. Finally, for large 
distances (r > z r ) the fifth dimension is manifestly com- 
pact, with an exponential cutoff in the correction to the 
potential. 

This general shape of the correction A(r) will be the 
same for all choices of the parameters /j, and z r . The only 
thing that will change when varying z r (and keeping /i 
fixed) is the position where the exponential cutoff occurs. 
Of course, if \xz r < 1 we will only see two different re- 
gions, as in this case the branes are so close together that 
the bulk space between them is essentially flat. 

log A 




FIG. 1: The figure shows the relative correction A to the 
Newtonian potential from l|27|l and for short and large 
distances, respectively, where V(r) = Vo(r)(l + A). The con- 
formal distance z r to the second brane is chosen such that 
fiz r — 10 15 . The figure clearly shows the transition from the 
one-brane result to that of a compact extra dimension around 
fir ~ 10 15 , i.e. around r ~ z r . 



V. SUMMARY 

In this paper we have considered the Randall-Sundrum 
model with two branes. The goal has been to study the 
effect on the Newtonian gravitational potential due to the 
finite distance between the branes. As expected, when 
the distance r in the potential is small compared to the 
brane separation z r , the potential is asymptotically the 
same as in the one-brane model where z r — > oo. In a 
sense, a short-range measurement of gravity is unable to 
see what is going on at a much larger scale. 

It is not until the distance r becomes comparable to z r 
that we can see the effect of the finite distance between 
the branes, and the potential changes from that of the 
one-brane model to that of a model with an ordinary 
compact extra dimension, where the relative correction 



to the four-dimensional potential decreases exponentially. 
The entire bulk-space dimension thus becomes invisible 
to a long-range experiment in gravity. 

Both of these examples illustrate a quite general fea- 
ture of almost all of physics - that a measurement done 
at a particular scale can only reveal information about 
the underlying physics at the approximately same scale. 
It should be noted, however, that it is not the physi- 
cal distance y r to the second brane that sets this scale, 
but rather the conformal distance z r ~ e Wr / /i, since all 
graviton masses are suppressed by the same factor e^ VT . 
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APPENDIX A: THE LIMIT fiz r < 1 

When the conformal distance to the second brane is 
small compared to the curvature radius of the fifth di- 
mension, we would expect the fifth dimension to behave 
like an ordinary compact dimension at all scales. More 
precisely, we should get the same potential V(r) for all 
distances r as with flat, five-dimensional space, where one 
of the dimensions is compact with an extension L — 2z r . 

When [iz r <C 1 the conformal distance z r between the 
branes is equal to the physical distance y r , since l+fj,z r — 
~ l+/j,y r . In this limit even the first mass eigenvalue 
will be very large, i.e. mi/jtt 3> 1, and we can use the 
asymptotic expressions for all Bessel functions in i|14[l . 
This yields tan(S - 2j£) ~ tan(^ - 2p + mZr ), with the 
solution 

m„ ~ — , [iz r <c 1 . (Al) 

This result should be compared to l|24|) . For large values 
of n the two expressions are identical to the lowest order, 
i.e. when l|24|l is expanded in powers of l/fiz r . However, 
in l|24|) it is better to use — ~ than m n ~ ^ even 

when fiz r 3> 1, because the next correction in the latter 
case would be ~ l//i 2 z 2 , whereas the next correction in 

(HP is ~ 1/(1 + pLZ r f < l/^ 2 z 2 . 

The quantization I|A1|) is exactly what we get from a 
flat, compact extra dimension, where the wavefunction is 
u{z) ~ e lmz and the periodic boundary condition gives 

e imL = lj with L = 2z r . 

With m n /fi > 1 we get |M(m,0)| 2 ~ l/z r just as in 
section llV Bl and the correction A to the Newtonian po- 
tential is therefore 

A ~ o ' 9 Y e- n ' r ^ ~ — -f , (A2) 

where we have used y r ~ z r . Again we get the same result 
as with a flat, extra dimension (see [lj). For short and 
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large distances, respectively, (|A21 is simplified to 

!2z r 
nr' (A3) 
2e-^/ z ", r»z r . 

Although a little more camouflaged than in the limit 
fiz r 3> 1, the result for r <C z r is once again the same as 
for a single brane, since the factor z r is absorbed into the 



five-dimensional Planck mass. Alternatively, if we had 
kept the factor (f - e~ 2 ^)/(l + ie~ 2 «^) in i|Al)l. we 
would get the exact same result as (12711 for r <C z r . 

In any case, the conclusion is that in the limit r <C z r 
the two-brane model is indistinguishable from the one- 
brane model, whereas in the limit r 3> z r the extra 
dimension behaves like an ordinary, flat-space compact 
dimension. 
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At this point, we are simply ignoring the radion field 
4>{x) ~ ha. The component h^, on the other hand, can 
be set to zero as a gauge choice 0,0. The radion will 
later be inserted back into the final result, using pd|. 
One can show that 

3 i 



h = uY2(u)J2(u)du — 
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For the definition of G, see e.g. the help system of Math- 
ematica (http://www.wolfram.com/ ). 



